Weak-field limit of /(i?)-gravity in three and more spatial dimensions 
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We investigate a point-like massive source in non- linear f(R) theories in the case of arbitrary 
number of spatial dimensions D > 3. If D > 3 then extra dimensions undergo toroidal compactifi- 
cation. We consider a weak-field approximation with Minkowski and de Sitter background solutions. 
In both these cases point-like massive sources demonstrate good agreement with experimental data 
only in the case of ordinary three-dimensional (D = 3) space. We generalize this result to the case 
of a perfect fluid with dust-like equations of state in the external and internal spaces. This perfect 
fluid is uniformly smeared over all extra dimensions and enclosed in a three-dimensional sphere. In 
ordinary three dimensional (D — 3) space, our formulas are useful for experimental constraints on 
parameters of f(R) models. 
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I. INTRODUCTION 

The idea of the multidimensionality of our Universe 
demanded by the theories of unification of the funda- 
mental interactions [l| is one of the most breathtaking 
ideas of theoretical physics. It is very important to sug- 
gest experiments which can reveal the extra dimensions. 
On the other hand, if we can show that existence of 
the extra dimensions is in contrast to observations, then 
these theories are prohibited. This important problem is 
extensively discussed in recent scientific literature (see, 
e -g-> [2]-[3)- in our paper we considered a point-like 
massive source in Kaluza-Klein models with an arbitrary 
number of toroidal internal spaces. The gravitational 
part of the action was taken in the linear (with respect 
to the scalar curvature R) form. It is well known that in 
general relativity in the weak-field limit, this is a good 
physical approximation to describe the famous gravita- 
tional experiments (perihelion shift, light deflection and 
time delay of radar echoes) . We expected that in the case 
of extra dimensions we can also satisfy the gravitational 
experiments if the sizes of the internal spaces will be small 
enough. To our surprise, it is not the case. We found 
that the point-like mass approach contradicts the obser- 
vational data if the total number of spatial dimensions 
D is more than three and this result does not depend on 
sizes of the internal spaces. Analysis performed in Refs. 
[HI [H| shows that the reason lies in equations of state. In 
the weak-field limit, point-like massive sources have dust- 
like equations of state in all spatial dimensions. How- 
ever, as it follows from our investigations, for considered 
Kaluza-Klein models, the concordance with observations 
requires non-zero equations of state in the internal spaces. 
It was shown in Ref. that the only compact astro- 
physical objects, that satisfy the observational data with 
the same accuracy as general relativity, are the latent 
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solitons with a dust-like equation of state in our three 
dimensions and definite (non-zero) equations of state in 
the internal spaces. Moreover, the condition of stabil- 
ity of the internal spaces singles out black strings/branes 
from the latent solitons and leads uniquely to p, t = — e/2 
equations of state in the internal spaces, and to the num- 
ber of the external dimensions do — 3. The main problem 
with the black strings/branes is to find a physically rea- 
sonable mechanism which can explain how the ordinary 
particles forming the astrophysical objects (e.g., Sun) can 
acquire rather specific equations of state (p* = —e/2) in 
the internal spaces. At the moment, it looks very prob- 
lematic. 

It is possible that a point-like massive source (with 
dust-like equation of state in external and internal 
spaces) does not contradict the observational data in 
some other generalized gravitational theories different 
from the one described in In this case we can 

circumvent the problem mentioned above. To clarify 
this question, in the present paper we consider non-linear 
f(R) theories of gravity for an arbitrary number of spatial 
dimensions D > 3. Starting from the pioneering paper 
[l2l |. the nonlinear theories of gravity f(R) have attracted 
the great deal of interest because these models can pro- 
vide a natural mechanism of the early inflation. Nonlin- 
ear models may arise either due to quantum fluctuations 
of matter fields including gravity [13], or as a result of 
compactification of extra spatial dimensions [l4| . Com- 
pared, e.g., to other higher-order gravity theories, f(R) 
theories are free of ghosts and Ostrogradski instabilities 
[IH . Recently, it was realized that these models can also 
explain the late-time acceleration of the Universe. This 
fact resulted in a new wave of papers devoted to this 
topic (see, e.g., reviews [Tr3|-[l9j). 

Therefore, in this article we investigate a point-like 
massive source in non-linear f{R) theories in the case 
of an arbitrary number of spatial dimensions D > 3. In 
a weak-field approximation these sources result in per- 
turbations of a background metrics. We consider two 
types of the background metrics: the Minkowski met- 
rics and the de Sitter metrics. Concerning the form of 
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the function f(R), we demand only that it should be 
an analytical function that can be expanded in a Tay- 
lor series in the neighborhood of the background solu- 
tions. The main result is that in both these cases point- 
like massive sources demonstrate good agreement with 
experimental data only in the case of ordinary three- 
dimensional (D — 3) space. We generalize this result 
to the case of a perfect fluid with dust-like equations of 
state in the external/our and internal spaces. This per- 
fect fluid is uniformly smeared over all extra dimensions 
and enclosed in a three-dimensional sphere. In ordinary 
three dimensional (D = 3) space, our formulas are use- 
ful for experimental constraints on parameters of f(R) 
models. 

The paper is structured as follows. In Section 2 
we consider a weak- field limit of f(R) theories with a 
point-like mass in the case of the Minkowski background 
and demonstrate good agreement with experimental data 
only for three-dimensional space. The similar result takes 
place in the case of the de Sitter background (see Sec- 
tion 3). In Appendix, we generalize our investigation to 
the case of perfect fluid enclosed in a three-dimensional 
sphere. The main results are summarized in the conclud- 
ing Section 4. 



one and we will investigate it in the weak- field limit. This 
means that the gravitational field is weak and the veloc- 
ities of test bodies are small compared to the speed of 
light c. In this case the metrics is only slightly perturbed 
from its flat spacetime value: 

9ik ~ Vik + h ik , (3) 

where are corrections of the order 1/c 2 . Then, up to 
the order 1/c 2 , the covariant components of Ricci tensor 
read [<|: 

R . k ^ 1 ( d 2 h\ | d 2 hj d 2 h\ d 2 h lk \ _ 

1 2 \dx k dx l dx l dx l dx l dx k dxWx 1 ) 

(4) 

We are going to investigate the weak-field approximation 
where the gravitational field is generated by a point-like 
mass at rest. To hold on the right-hand side of Eqs. ([1]) 
and ([2]) the terms up to the order 1/c 2 , components of 
the energy-momentum tensor are approximated as 

T 00 «pc 2 , T 0Q «0, T q/3 «0 =► T = T ik g ik » pc 2 , 
where the rest mass density is 



II. WEAK-FIELD LIMIT OF /(ii)-GRAVITY. 
FLAT SPACETIME BACKGROUND 

It is well known (see, e.g., [2(| HH) that in the case 
of f(R) gravitational theories with an arbitrary number 
of spacetime dimensions V = 1 + D > 4, the Einstein 
equations read 

f'{R)Rik - \f{R)gik - [f'(R)h,k 



-9ik{f'(R)U;n9 r ' 



■T, 



ik •■ 



(i) 



where S D = 2n D / 2 /T(D /2) is the total solid angle (the 
surface area of the (D— l)-dimensional sphere of a unit ra- 
dius), Gv is the gravitational constant in D-dimensional 
spacetime, the prime denotes differentiation with respect 
to the scalar curvature R = Rik9 %k '- f'(R) = df/dR, 
and the semicolon ; denotes the covariant derivative with 
respect to the metric coefficients g^. The trace of Eq. 
© is 



f'(R)R 



l + D 



f(R)+D[f(R)}, m , n g^ 



2SdGx 



■T . 



(2) 

Hereafter, the Latin indices i,k = 0,...,D, the Greek 
indices a, j3 = 1, . . . , D and we use the sign convention for 
the metrics, the Riemann and Ricci tensors in accordance 
with the book (22 1. 

In this section, we assume that in the case of absence 
of a matter source the spacetime is the Minkowski space- 
time: goo = %o = 1, 5oq = VQa = 0, g a p = r) a p = —S a p. 
In presence of matter, the metrics is not a Minkowskian 



p = mS(r) 



(6) 



and r = (a; 1 , x 2 , . . . , x D ) is a D-dimensional radius vec- 
tor. 

Now, let us suppose that f(R) is an analytical func- 
tion and can be expanded in a Taylor series near the flat 
spacetime background value R = : 

f(R) = R + aR 2 + o(R 2 ) , \R\ > \a\R 2 , (7) 

where o (R n ) means that this function contains the terms 
R m with powers m > n. In Eq. ©, a = (1/2) f" (0) 
and we assumed that /(0) = 0. The latter means that 
the cosmological constant A is absent in the model in 
accordance with our assumption that the background 
spacetime is flat. We also normalize our model in such 
a way that /'(0) = 1 which provides the transition to 
the usual form of the Einstein equations in the limit 
R -> f(R) -> R. From Q we get 



f'(R) = l + 2aR + o(R) 



(8) 



Substituting and © into ([T]), we obtain (up to the 
order 1/c 2 ) 



Rik — 2aR-i-k ~ — D , V Tik + -Rrjik — ^arjikR-.m-.nV 7 ' 
c 1 2 

Taking into account that up to the order 1/c 2 

d 2 R 



R-.i-.k 



dx l dx k 



(9) 
(10) 
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we can write Eq. ([9]) in the following form: 

i ( d 2 h[ . d 2 h{ d 2 h\_ _ iX d 2 h lk 



\ dx k di 



■ l ' dx l dx l dx l dx k ^ dxWx 1 



-2a Vik ~™» 



dx m dx n 



V 



(11) 



where we used Eq. (jj). With the help of the gauge 
condition 



d 

dx k 



(h k - l -(h\+±aR)8^^Q (12) 
the formula (JTTJ) can be written in the form 



\ Ah lk « ——j-—Tik + ^-Rijik + 2arj lk AR , (13) 
2 cr 2 



where A = 5 a/3 - 



i:) 2 



is the D-dimensional Laplace 



dx a dxP 

operator and up to the order 1/c 2 holds 

We now turn to Eq. $2$. Up to the order 1/c 2 , Eq. (J2j) 
reads 

D — 1 . 2S nGjy 
R-2aDARm , V T. 15) 

To get it, we used Eqs. Q, ©, (dUJ) and CG!}. Combining 
Eqs. (JT3J) and (|T5]). we obtain 

1 a , 25 D Gx, / 1 \ 2a 

-Aftifc « — ^ — \ Tik - D1 Tr) ik I - - _ ^ kAR. 

(16) 

It is worth noting that the weak- field equations (|13l) , ([15]) 
and (I16[) are valid for any energy-momentum tensor Ti k ~ 
0(c 2 ). Eq. (fT6|) shows that for a dust-like perfect fluid 
(|5|) (with an arbitrary form of the mass density p = p(r)) 
we have 

2(j) Aid 



/8 



where the function 1 ^ satisfies the D-dimensional Poisson 
equation 



A<j> = SdGvP 



(18) 



1 From the expression for /loo in (1171) , it is clear that the function <j> 
defines only a part of the non-relativistic gravitational potential. 
Another part follows from the i?-term. It originates from the 
additional scalar degree of freedom of non-linear f(R) theories. 



if the D-dimensional gravitational constants Gx> and G-p 
are related as follows: 



(19) 



Gt> — 2Gj)— — - 



The solutions of the Z?-dimensional Poisson equation (|T8|) 
in the case of a point-like mass source were obtained in 
[lil Hi]. To solve (HU) we should specify the topology of 
space and the boundary conditions. As for the boundary 
conditions, we require that the non-relativistic gravita- 
tional potential (p(r) should have the newtonian behav- 
ior ~ 1 /ra , where is the length of a radius vector in 
three-dimensional space, at far distances from the gravi- 
tating body. To be more precise, we require the following 
boundary condition: 



lim (r 3 h 00 ) 

T-3-S- + 00 



•9 ' 



(20) 



where r„ 



- 2Gatto/c 2 is the three-dimensional 
Schwarzschild radius of the gravitating body with the 
mass m and Gn is the Newtonian gravitational con- 
stant. In the case of non-compact extra dimensions, i.e. 
when the (D = 3 + d)-dimcnsional space has topology 
Mn = R 3+rf , the function behaves as ~ l/r 1+d , where 
r is the length of a radius vector in D-dimensional space 
[23l |24T ] . Obviously, it contradicts the boundary condi- 
tion ((201) ■ Therefore, we suppose that the (D = 3 + <£)- 
dimensional space has the factorizable geometry of a 
product manifold M D = M 3 x T d . Here R 3 describes 
the three-dimensional flat external (our) space and T d 
is a torus which corresponds to a d-dimensional internal 
space with volume Vd — Y\.i=x a ii where a, are the periods 
of tori. For this topology and with the boundary condi- 
tion that at infinitely large distances from the gravitating 
body the potential must go to the Newtonian expression, 
we can find the exact solution of the Poisson equation 
(|T5]) [23l [24} . The boundary condition requires that the 
multidimensional Gt> and Newtonian Gjy gravitational 
constants are connected by the following condition: 



S D G V /V d = 4ttG 



jv ■ 



(21) 



Assuming that we consider the gravitational field of the 
gravitating mass m at distances much greater than peri- 
ods of tori, we can restrict ourselves to the zero Kaluza- 
Klein mode. For example, this approximation is very well 
satisfied for the planets of the solar system because the 
inverse-square law experiments show that the extra di- 
mensions in Kaluza-Klein models should not exceed sub- 
millimeter scales (25|. Then, the solution of Eq. (fT8|) 
reads 



(r) 



Gntti 



T3 



2r 3 



(22) 



In the case when the gravitating mass is uniformly 
smeared over the extra dimensions, the approximate for- 
mula (|2"2"]) becomes the exact equality [23|, [2J] . 
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Let us turn now to Eq. (|15j) for the scalar curvature 
R. For a point-like mass, this equation takes the form 



AR- h 2 Rk. - 



aDc 2 



<J(r) ee -P6(r) , (23) 



where 



D-l 

AaD 



1/2 



>0. 



(24) 



Obviously, /i plays the role of a mass for the additional 
scalar degree of freedom in /(i?)-gravities. To have a 
physically reasonable (stable) solutions, we should de- 
mand > 0. It happens for a = (1/2)/" (0) < 0. Some 
papers require the opposite sign for f"(Ro). However, it 
can be easily realized that the choice of this sign depends 
on the sign convention for the metrics and the definition 
of the Riemann tensor. The authors in Refs. [26l - [28| 
have chosen the sign convention similar to ours and also 
found that f"(Ro = 0) should be negative. In (D = 3)- 
dimensional space we get 



1 



3|/"(0)| ' 



(25) 



which coincides with the well known expression for the 
mass of the scalar degree of freedom in the case of 
the background solution Rq = and the normalization 
/'(0) = 1 (see, e.g., [HE3). 

We consider first the case when space has topology 
Md = M 3+d , i.e. the extra dimensions are non-compact. 
Then, the decreasing solution (i.e. lim R = 0) of Eq. 

T— >+oo 

(|23|) in the region r > is 



1/2 



1 — D 

r 2 e" 



R^CKo^ry 1 - 

where K^ D / 2 )-i is the modified Bessel function of the 
second kind and C is the constant of integration. In par- 



-1/2 



ticular case D = 3, we get R ~ CK 1 / 2 (fj,r 3 )r 3 

1 /2 

C( 7 r/(2/i)) i/ e-^/r-s that gives the Yukawa contribu- 
tion to the Newtonian gravitational potential. It can 
be easily seen that in the case of non-compact extra 
dimensions {D > 3), neither function 4> nor function 
R provide the boundary condition (|2T)|) . Therefore, it 
is necessary to consider space with topology Mn — 
R 3 x T d and compact extra dimensions. For this topol- 
ogy, it is natural to impose periodic boundary condi- 
tions: fl(r 3 , £i,£ 2 , •-•,&) = -Rfo.fi, + 
Oi, . . . , £d), i = 1) • • • , d, where ai denotes a period in the 
direction of the extra dimension £^ . Then it is convenient 
to rewrite Eq. (1231) as follows: 

a * r + e - ^ 2r « -^( r 3) n 



1 +2? (2nk\ 

J fc=-oo v ' 



(27) 



We seek a solution in the form 

<l ^ +00 "x. 

R " 



d ^ +00 +00 

Ilr E • E Ai...*„( r 3) 

3 h,— - 



3=1 



'27rfci 
X cos I £1 



00 kd — — oo 

( "Znkd 

. cos 



where the function /fc 1 ...fe d (r 3 ) satisfies the equation 

2 



(28) 



Az]hx...k d 



d 

E 

i=i 



27rfc 1 



fk!...k d = -P6(r 3 ). 



The only solution of this equation with the boundary 





condition lim /fe 1 ...fc d (r 3 ) = is 

r 3 ->+oo 



fk t ...k d (r 3 ) 



x exp 



47TT3 




Taking into account that for r 3 
kj = 0, j = 1, d gives the main contribution, we obtain 
at large distances the following asymptote 



R 



D - 1 



3=1 



a, 47rr3 



AaD(D - 2) 



?"3 



where we took into account the relation (|2Tj) between the 
Newtonian Gat and multidimensional Gx> gravitational 
constants. 

It is worth noting that in the case of the smeared extra 
dimensions (i.e. when the gravitating mass is uniformly 
smeared over all extra dimensions), Eq. ([27| is reduced 
to 

d , 

A 3 R-^R^ _/3(J(r 3 )TJ — , (32) 

3=1 a * 

and the only solution of this equation with the boundary 
condition lim R = is the function (f3"T|) . 

Thus, the substitution of expressions (|22|) and (|3T)) in 
Eq. (|17|) results in the following form of hoo and h a p- 



•_9_ 



1 



D(D-2) 
1 1 



-fJ.r 3 



rz \D-2 



(33) 

e~^ 3 ) S a p . (34) 



D(L> - 2) 

With the help of the Green function similar expressions 
were obtained in (2|| in the particular case D = 3. As 
it follows from (|33[) , the non-relativistic gravitational po- 
tential is 



c 2 G N m 
¥V3 = -^-n-oo = 

2 ?3 



1 



1 



D(D -2) 



A* 



= |3/"(0)| 1/2 > 0. 



(35) 



15=3 
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Thus, the additional scalar degree of freedom results in 
the Yukawa correction to the non-relativistic gravita- 
tional potential. Neither the inverse square law (ISL) 
experiments nor the solar system gravitational tests re- 
veal so far such corrections. They give only upper limits 
on the parameters of the Yukawa interaction. For exam- 
ple, in the case D = 3 in the formula the Yukawa 
parameter a (which is the prefactor in front of the expo- 
nent) is equal to a = 1/{D(D - 2)) = 1/3 « 1CT 1 / 2 that 
results in the upper limit [25j 

A < A m ax ~ 7xlCT 3 cm => |/"(0)| < 1.6x l(T 5 cm 2 . 

(36) 

A similar restriction on /"(0) was also obtained, e.g., in 
[28l l3l| . In the terms of the mass of the scalar degree of 
freedom it gives /i > 2 x 1CP 3 eV. These estimates show 
that r 3 ^> A for the classical gravitational tests (e.g., de- 
flection of light) in solar system, and we can neglect the 
Yukawa corrections at these distances. Therefore, for 
D > 3, we arrive at a discrepancy between hoo ~ —r g /r^ 
and h aa ~ —(1/(1? — 2))r g /r^,. In linear (with respect to 
R) Kaluza-Klein models with a point-like massive source, 
exactly this discrepancy leads to a contradiction with the 
classical gravitational tests |9j . Therefore, for considered 
in this section non-linear f(R) models 2 , a point-like mas- 
sive source demonstrates good agreement with experimen- 
tal data only in the case of ordinary three-dimensional 
(D = 3) space. In Appendix we demonstrate that the 
similar conclusion holds also in the case of a dust-like 
spherically symmetric compact gravitating source. 



III. DE SITTER SPACETIME BACKGROUND 

According to current observational data, dark energy 
dominates the energ y d ensity in the modern Universe 
(see, e.g., the lecture [32|])- There is a number of theoret- 
ical approaches that have been adopted so far to explain 
the origin of dark energy [33j . It is usually assumed that 
the Universe at late times enters the de Sitter stage. It 
is well known that the de Sitter space is the constant 
curvature space: 



Roi. 



D - 1 



A.g ifc 



Ro 



2{D + 1) 
1 



A. 



(37) 

where A is the cosmological constant, and according to 
the observations A is extremely small at present time: 
A ~ 10~ 57 cm -2 . Consequently, in contrast to the previ- 
ous section, we suppose that the background space is de 
Sitter space, rather than Minkowskian one, with a back- 
ground metrics goik and the background curvature Rq. 
Dealing with f(R) theories, we also assume that Ro de- 
fines a stable de Sitter fixed point of f(R) theory and A is 



connected with f(Ro) (see Eq. (|43|) below). The formula 
(f3"T|) shows that Ro is the same order of magnitude as A 
and should be much less than the curvature R = Ro + Ri 
in the vicinity of astrophysical objects, i.e., \Ri\ 3> |-Ro|- 
Note that we still consider the weak-field approximation. 
Then, instead of the expansion ([7]) we have 

f{R) = f{R ) + f'(R )Ri + \f"{Ro)R\ + o{R\) . (38) 

We also suppose for expansions f{R) and f'{R) the fol- 
lowing inequalities: 



f{Ro) + f(R )R 1 » -jW(Ro)I%, n = 2, 3, ... , 
n ! 

f'(Ro) + f"(R )Ri » ^/ (n+1) (^o)i?r, n = 2, 3, ... . 

(39) 



Taking into account Eqs. (|38| and (|39|) . we obtain from 
Eq. © 

f'{Ro)Roik + f"(Ro)RlRoik + f'{Ro)Rlik 
-\f{Ra)9ik - \f'{.Ro)Rl9ik - [f'(Ro)];i;k 



-[f"(R )Rihk + g*k[f'(Ro)U;ng r ' 

+9ik[f"(Ro)Rl]; m ; n 9" 



k-nTj 



V-Lik 7 



(40) 



where R\,Ruk are the differences between total per- 
turbed values R,Rik and corresponding background val- 
ues Ro,Roik, and we introduced the notation k-p = 
2SdGv I 'c 4 . Here, we do not require that R\, Ruk are 
defined up to the linear terms guk — gik — goik ■ If nec- 
essary, the corresponding linearized expressions can be 
found, e.g., in (34J. In what follows, we consider the case 
when background value R = const. Without loss of 
generality, we also suppose that f'(Ro) = 1- Under these 
conditions, we get from (|4"U]) 

Roik + f" {Ro)RlRoik + Rlik — -^f{Ro)goik 

--^fiR^giik - -^Rigoik - f" ' {Ro)[Ri];i-k 



+goikf"(Ro)[Ri] 

Zero-order approximation gives 
Roik - -^f(Ro)goik = = 



kvTik ■ 



(41) 



Ro — 



l + D 



f(Ro) , 
(42) 

which demonstrates that f(Ro) plays the role of the cos- 
mological constant (cf. with (|37l0 : 



/(-Ro) 



D - 1 



-A. 



(43) 



We require only that the function f(R) can be expanded in series 
0. 



Obviously, the parameter r = —Rof'(Ro)/f{Ra) (de- 
fined in J35|]) is equal to tq = — 2 for D — 3 what should 



G 



hold for the de Sitter fixed point. Taking into considera- 
tion Eqs. (f4*2"j). we rewrite Eq. (|4*Tj) as follows 



f"(R )RiR 



1 + D 



-goik 



Rq 

l + D 



giik + Rlik 



--^Rl90lk — f" ' (Ro)[Rl];i;k 

-goikf" {Ro)[Ri] 



(44) 



Evidently, because of extreme smallness of A, the back- 
ground metrics can be approximated by the Minkowskian 
one: goik ~ Vik 3 - For the same reason, in solar system 
(e.g., in vicinity of Sun) \Rq\ <C Therefore, we can 

drop two first terms, containing i? , in the left hand side 
of Eq. (|44l) and replace goik by r\i k . This reduces Eq. 
(|4"4"| to Eq. ©. The only difference is that a is equal 
now to f"{R )/2. 

Let us consider now Eq. @. In this case, in complete 
analogy with the derivation of Eq. (l40l) , we obtain 

f'{Ro)R + f"(Ro)RoRi + f'(Ro)Ri - ^^-f(Ro) 
-^^f(Ro)Ri + D[f'(R )], m , n g mn 



-D[f"(Ro)R 1 ], m , n g mn ^ 2 -^^T . 



(45) 



It can be easily seen that in zero-order approximation 
we again obtain Eq. (|42p . Taking into account (l4"2"1) . the 
normalization f'(Ro) = 1, the constancy of Ro and the 
decomposition g ik = g Qlk + g lik , we get 



[Rl]-m;ng° mn 

~ c±Df"{R a 



D - 1 
2Df"(R ) 



D Ro 



Ri 



■T. 



(46) 



Therefore, the mass squared of the scalar degree of free- 
dom is 



D - 1 
2Df"(R ) 



15 R ° 



(47) 



In the case of three-dimensional space D — 3 and the 
condition f'(Ro) = 1, the similar formula for the mass 
squared can be found, e.g., in [2^, |30|, HI]. In the ap- 
proximation |i? | *C we can drop in Eqs. (|45l) - (|47|) 
the terms proportional to Rq. Then, keeping in mind the 
relation goik ~ Vik, we can easily see that Eq. (|46|) is 
reduced to Eq. (fT5|) where /"(0) should be replaced by 
/"(i?o). Therefore, up to this replacement, the formulas 
(f3"3"| and (|3~4"|) and following from them conclusions are 



3 For example, in the case D = 3 this approximation works very 
well at distances r3 <C A -1 / 2 ~ 10 28 cm (see, e.g., the formula 



also suitable in the case of the de Sitter background 4 . 
For example, the ISL experiments require: 



\f"(Ro)\ < 1-6 x l(T 5 cm 2 



(48) 



This estimate indicates that {/"(Rq)^ 1 > 10 5 cm~ 2 > 
i?o| ~ A ~ 10 _57 cm -2 and we can really drop the terms 
proportional to R in Eqs. P5|) - (|47)l . 

Therefore, in this section we arrive at the same con- 
clusion as in the previous section: for non- linear f(R) 
models satisfying the conditions (j3"5)) . (|39l) and with the 
de Sitter background, a point-like massive source demon- 
strates good agreement with experimental data only in 
the case of ordinary three-dimensional (D — 3) space. 



IV. CONCLUSION 

It is well known that a point-like mass is a very good 
approximation for classical gravitational tests in general 
relativity with three spatial dimensions. In our paper 
we checked this approximation for non-linear f(R) theo- 
ries for an arbitrary number of spatial dimensions D > 3. 
We performed our investigations in a weak-field limit and 
considered two forms of the background solutions when 
a gravitating matter is absent. First, it is the Minkowski 
metrics. Second, keeping in mind dark energy in the Uni- 
verse, it is the de Sitter metrics. Concerning the form of 
the function /(i?), we demand only that it should be 
an analytical function that can be expanded in a Taylor 
series in the neighborhood of the background solutions. 
For both of these cases we obtained the similar expres- 
sions for the perturbations of the metric coefficients hoo 
and h a /3- This makes it possible to calculate, for exam- 
ple, the deflection of light in solar system. The form 
of the perturbations hoo and h a p exactly indicate that 
for both background metrics we get the same conclusion: 
point-like massive sources demonstrate good agreement 
with experimental data only in the case of ordinary three- 
dimensional (D = 3) space. We generalize this result to 
the case of a perfect fluid with dust-like equations of state 
in the external/our and internal spaces. This perfect fluid 
is uniformly smeared over all extra dimensions and en- 
closed in a three-dimensional sphere. It is not difficult to 
see that more complicated form of the internal structure 
does not change the situation. It follows from the fact 
that the asymptotic ratio h aa /hoo = l/(-D — 2) (which 
is in agreement with the observations only for D = 3) 
takes place for an arbitrary form of p{r) (see Eq. (jTTJ)) . 



(5.76) in [T|). 



4 It is worth noting that despite the fact that in our approach 
the equations with the de Sitter background look similar to the 
equations with the flat background, the replacement of /"(0) 
by f"(Ro) allows us to investigate a class of models with /" 
divergent at i?o = 0. For example, in models of the form of 
f(R) = R 1+c /[(l + e)R e ] , e < 1 (see, e.g., Ref. [H) we get 
f"(Ro) = e/R<3, and f"(Ro) can be made arbitrary small for 
any small non- zero values of _Ro- 
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We think that the found contradiction with the obser- 
vations is a generic property of the Kaluza-Klein mod- 
els. The point is that in the case of dust-like equations 
of state in the external and internal spaces, both point- 
like and extended gravitating objects result in the men- 
tioned above asymptotic ratio between the metric per- 
turbations: haa/hoo = 1/(D — 2). It happens, e.g., in 
Kaluza-Klein models with toroidal 0-[Hl and spherical 
38] compactifications as well as in considered above f(R) 
models. It is also clear that the exterior to a black hole 
may not coincide with the field exterior to a point mass. 
We can see this on the example of the soliton solutions 
in Kaluza-Klein models with toroidal compactification. 
Here, latent solitons, black strings and black branes do 
not contradict the observations (see, e.g., [13, EH)- How- 
ever, in all these cases we get very strange equations of 
state in the internal space. It is of interest to obtain an 
analogue of the black string/brane solutions for multidi- 
mensional f(R) models to see what the form of equations 
of state takes place here. We are going to undertake such 
a study in the near future. 

In conclusion, we would like to note that in ordinary 
three dimensional D = 3 space, the formulas obtained 
in our paper are useful for experimental constraints on 
parameters of f(R) models. For example, we found the 
formula for the gravitational force between two spheres. 
This formula can be used for calculation f"(Ro) from 
direct experimental measurements. 
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constants. These equations have the following solutions: 
Gniti 



7*3 

2-KG N pzrl 



2irG N p 3 r s 



2rl 



2r.„ 



(A2) 



To get these solutions, we use the boundary condi- 
tions 4>out( r 3 — ► +oo) — > 0, |0m(O)| < +oo, and the 
matching conditions </w(r s ) = 4> ln {r s ), d(f> out /dr 3 \ ra = 
d<j> in /dr 3 \ rs . 

Similarly, Eq. (|T5j) outside and inside the sphere takes 
the form 



A3 Rout — 
A 3 i?i„ 



P?Rout 



V R% 



1 f 

a 

1 d 2 . 
— -T2\r?,Ri 
r 3 dr% 



T^yTaRout) 

r 3 dr% 



P- 2 Rout 



0. 



P?R V 



where the mass squared p 2 is defined by Eq. 
introduce the notation 



P~3 = Pp3 , P 



1 



D - 1 



2aD(D - 2) 



AnG N p 3 , 
(A3) 
and we 

(A4) 



Using the boundary and matching conditions similar to 
conditions for the functions <f>, we obtain the following 
solutions: 

e -A«-3 4TrG N p 3 
Rout = 3 |smh(/zr s ) - pr s cosh(/xr s )] 



T3 

-Gmiti 



r 3 



/i 3 rf 



[pr s cosh(pr s ) - smh(jj,r s )} 



(A5) 



V. APPENDIX: DUST-LIKE PERFECT FLUID 
IN A THREE-DIMENSIONAL SPHERE 

In this appendix we consider the case where a point- 
like mass is replaced by a dust-like perfect fluid with 
the energy- momentum tensor ([5]). We suppose that 
this perfect fluid is uniformly smeared over all extra di- 
mensions and, at the same time, enclosed in a three- 
dimensional sphere of the radius r s . Then, multidimen- 
sional p and three-dimensional p 3 energy densities are 
related as p 3 — pVd- We also assume that p 3 = const. 
Therefore, Eq. (fT8|) outside and inside the sphere, re- 
spectively, reads 



A 



3f>out 



A 



1 d 2 

-JJ^^out) = 0, 
1 d 2 



3Wn 



-T2( r 3( 

r 3 dr% 



4itGnP3 



(Al) 



where we take into account the relation (|2"Tj) between 
multidimensional Gt> and Newtonian Gn gravitational 



and 



Rir>. 



4iTG N p 3 
p 2 



e ' irs (1 + pr s ) sinh(^r 3 ) - 1 

pr 3 



With the help of the asymptotic formula a; cosh a; — 
sinhx w x 3 /3 for x <C 1, we can easily see that R ou t 
is reduced to Eq. (f3~Tj) for r s — > 0. On the other hand, 
up to the prefactor (3, Ri n and R ou t are reduced to 4n n 
and (f>out (|A2|) . respectively, for p, — > 0. Now, we can get 
the metric coefficients /loo and h a p defined in Eq. (|17p . 
For example, we obtain outside of the sphere 



r^out 

"■00 



LA 

r 3 



1 



1 



D(D - 2) 



x-3-o-[/w« cosh(^r s ) - sinh(/ir s )] 

,out _ _lg_ ( 1 1 

ap r 3 \D-2 D(D~2) 



(A7) 



x— — [pr s cosh(/ir s ) - sinh(^r s )] I S a/3 . (A8) 



p 3 r 3 
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Similar to a point-like massive source, we also arrive at 
the relation h%$ h°f/{D - 2) for ^r 3 > 1. 

It is worth noting that the non-relativistic potential 
Pout {r 3) = (c 2 /2)/inn* coincides formally with the for- 
mula (3.7) in [24| where we should keep the first Kaluza- 
Klein modes k = ±1 and take into account that /i = 
2ttx±1i R2 = f S) Ri = 0. We also need to replace in 
(3.7) the Yukawa parameter a = 2 (for simplicity we 
consider in (3.7) the case of one extra dimension) by 
a = 1/(D(D — 2)). We can apply this analogy to an- 
other formulas in [23 |. For example, the gravitational 
force between two spheres (with masses mi and m 2 and 
radii r s \ and r S 2) is described by the formula similar to 
the expression (4.7) in [24| : 



where F/v(r 3 ) = — GjvWima/rf. Obviously, it makes no 
sense to use this formula for D > 3 because this case, as 
shown above, contradicts the observational data. How- 
ever, we can apply it to the case D = 3 where a = 1/3. 

1 /2 

Therefore, since A| D=3 = |3/"(i?o)| , then we can 
use this formula for direct experimental measurement of 

f'(Ro)- 



F(r 3 ) 
F N (r 3 ) 



X 



= 1 + l a ( AV ( AV I* e -(^.i-r-)/A 
4 \r sl J \r s2 J A 



(A9) 
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